Traceability of SARS-CoV-2 transmission through quasispecies analysis.

During COVID-19 pandemic, consensus genomic sequences were used for rapidly monitor the spread of the virus worldwide. However, less attention was paid to intrahost genetic diversity. In fact, in the infected host, SARS-CoV-2 consists in an ensemble of replicating and closely related viral variants so-called quasispecies. Here we show that intrahost single nucleotide variants (iSNVs) represent a target for contact tracing analysis. Our data indicate that in the acute phase of infection, in highly likely transmission links, the number of viral particles transmitted from one host to another (bottleneck size) is large enough to propagate iSNVs among individuals. Furthermore, we demonstrate that, during SARS-CoV-2 outbreaks when the consensus sequences are identical, it is possible to reconstruct the transmission chains by genomic investigations of iSNVs. Specifically, we found that it is possible to identify transmission chains by limiting the analysis of iSNVs to only three well-conserved genes, namely nsp2, ORF3, and ORF7.

TaME-seq2: tagmentation-assisted multiplex PCR enrichment sequencing for viral genomic profiling.

BACKGROUND: Previously developed TaME-seq method for deep sequencing of HPV, allowed simultaneous identification of the human papillomavirus (HPV) DNA consensus sequence, low-frequency variable sites, and chromosomal integration events. The method has been successfully validated and applied to the study of five carcinogenic high-risk (HR) HPV types (HPV16, 18, 31, 33, and 45). Here, we present TaME-seq2 with an updated laboratory workflow and bioinformatics pipeline. The HR-HPV type repertoire was expanded with HPV51, 52, and 59. As a proof-of-concept, TaME-seq2 was applied on SARS-CoV-2 positive samples showing the method's flexibility to a broader range of viruses, both DNA and RNA. RESULTS: Compared to TaME-seq version 1, the bioinformatics pipeline of TaME-seq2 is approximately 40× faster. In total, 23 HPV-positive samples and seven SARS-CoV-2 clinical samples passed the threshold of 300× mean depth and were submitted to further analysis. The mean number of variable sites per 1 kb was ~ 1.5× higher in SARS-CoV-2 than in HPV-positive samples. Reproducibility and repeatability of the method were tested on a subset of samples. A viral integration breakpoint followed by a partial genomic deletion was found in within-run replicates of HPV59-positive sample. Identified viral consensus sequence in two separate runs was > 99.9% identical between replicates, differing by a couple of nucleotides identified in only one of the replicates. Conversely, the number of identical minor nucleotide variants (MNVs) differed greatly between replicates, probably caused by PCR-introduced bias. The total number of detected MNVs, calculated gene variability and mutational signature analysis, were unaffected by the sequencing run. CONCLUSION: TaME-seq2 proved well suited for consensus sequence identification, and the detection of low-frequency viral genome variation and viral-chromosomal integrations. The repertoire of TaME-seq2 now encompasses seven HR-HPV types. Our goal is to further include all HR-HPV types in the TaME-seq2 repertoire. Moreover, with a minor modification of previously developed primers, the same method was successfully applied for the analysis of SARS-CoV-2 positive samples, implying the ease of adapting TaME-seq2 to other viruses.

Temporal dynamics of SARS-CoV-2 genome mutations that occurred in vivo on an aircraft.

We analyzed variations in the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome during a flight-related cluster outbreak of coronavirus disease 2019 (COVID-19) in Shenzhen, China, to explore the characteristics of SARS-CoV-2 transmission and intra-host single nucleotide variations (iSNVs) in a confined space. Thirty-three patients with COVID-19 were sampled, and 14 were resampled 3-31 days later. All 47 nasopharyngeal swabs were deep sequenced. iSNVs and similarities in the consensus genome sequence were analyzed. Three SARS-CoV-2 variants of concern, Delta (n=31), Beta (n=1), and C.1.2 (n=1), were detected among the 33 patients. The viral genome sequences from 30 Delta-positive patients had similar SNVs; 14 of these patients provided two successive samples. Overall, the 47 sequenced genomes contained 164 iSNVs. Of the 14 paired (successive) samples, the second samples (T2) contained more iSNVs (median: 3; 95% confidence interval [95%CI]: 2.77-10.22) than did the first samples (T1; median: 2; 95%CI: 1.63-3.74; Wilcoxon test, P=0.021). 38 iSNVs were detected in T1 samples, and only seven were also detectable in T2 samples. Notably, T2 samples from two of the 14 paired samples had additional mutations than the T1 samples. The iSNVs of the SARS-CoV-2 genome exhibited rapid dynamic changes during a flight-related cluster outbreak event. Intra-host diversity increased gradually with time, and new site mutations occurred in vivo without a population transmission bottleneck. Therefore, we could not determine the generational relationship from the mutation site changes alone.

Accelerated SARS-CoV-2 intrahost evolution leading to distinct genotypes during chronic infection.

The chronic infection hypothesis for novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant emergence is increasingly gaining credence following the appearance of Omicron. Here, we investigate intrahost evolution and genetic diversity of lineage B.1.517 during a SARS-CoV-2 chronic infection lasting for 471 days (and still ongoing) with consistently recovered infectious virus and high viral genome copies. During the infection, we find an accelerated virus evolutionary rate translating to 35 nucleotide substitutions per year, approximately 2-fold higher than the global SARS-CoV-2 evolutionary rate. This intrahost evolution results in the emergence and persistence of at least three genetically distinct genotypes, suggesting the establishment of spatially structured viral populations continually reseeding different genotypes into the nasopharynx. Finally, we track the temporal dynamics of genetic diversity to identify advantageous mutations and highlight hallmark changes for chronic infection. Our findings demonstrate that untreated chronic infections accelerate SARS-CoV-2 evolution, providing an opportunity for the emergence of genetically divergent variants.

Viral Population Heterogeneity and Fluctuating Mutational Pattern during a Persistent SARS-CoV-2 Infection in an Immunocompromised Patient.

Literature offers plenty of cases of immunocompromised patients, who develop chronic and severe SARS-CoV-2 infections. The aim of this study is to provide further insight into SARS-CoV-2 evolutionary dynamic taking into exam a subject suffering from follicular lymphoma, who developed a persistent infection for over 7 months. Eight nasopharyngeal swabs were obtained, and were analyses by qRT-PCR for diagnostic purposes. All of them were considered eligible (Ct < 30) for NGS sequencing. Sequence analysis showed that all sequences matched the B.1.617.2 AY.122 lineage, but they differed by few mutations identifying three genetically similar subpopulations, which evolved during the course of infection, demonstrating that prolonged replication is paralleled with intra-host virus evolution. These evidences support the hypothesis that SARS-CoV-2 adaptive capacities are able to shape a heterogeneous viral population in the context of immunocompromised patients. Spill-over of viral variants with enhanced transmissibility or immune escape capacities from these subjects is plausible.

Assessing the hidden diversity underlying consensus sequences of SARS-CoV-2 using VICOS, a novel bioinformatic pipeline for identification of mixed viral populations.

INTRODUCTION: Coinfection with two SARS-CoV-2 viruses is still a very understudied phenomenon. Although next generation sequencing methods are very sensitive to detect heterogeneous viral populations in a sample, there is no standardized method for their characterization, so their clinical and epidemiological importance is unknown. MATERIAL AND METHODS: We developed VICOS (Viral COinfection Surveillance), a new bioinformatic algorithm for variant calling, filtering and statistical analysis to identify samples suspected of being mixed SARS-CoV-2 populations from a large dataset in the framework of a community genomic surveillance. VICOS was used to detect SARS-CoV-2 coinfections in a dataset of 1,097 complete genomes collected between March 2020 and August 2021 in Argentina. RESULTS: We detected 23 cases (2%) of SARS-CoV-2 coinfections. Detailed study of VICOS's results together with additional phylogenetic analysis revealed 3 cases of coinfections by two viruses of the same lineage, 2 cases by viruses of different genetic lineages, 13 were compatible with both coinfection and intra-host evolution, and 5 cases were likely a product of laboratory contamination. DISCUSSION: Intra-sample viral diversity provides important information to understand the transmission dynamics of SARS-CoV-2. Advanced bioinformatics tools, such as VICOS, are a necessary resource to help unveil the hidden diversity of SARS-CoV-2.

SARS-CoV-2 intra-host single-nucleotide variants associated with disease severity.

Variants of severe acute respiratory syndrome coronavirus 2 frequently arise within infected individuals. Here, we explored the level and pattern of intra-host viral diversity in association with disease severity. Then, we analyzed information underlying these nucleotide changes to infer the impetus including mutational signatures and immune selection from neutralizing antibody or T-cell recognition. From 23 January to 31 March 2020, a set of cross-sectional samples were collected from individuals with homogeneous founder virus regardless of disease severity. Intra-host single-nucleotide variants (iSNVs) were enumerated using deep sequencing. Human leukocyte antigen (HLA) alleles were genotyped by Sanger sequencing. Medical records were collected and reviewed by attending physicians. A total of 836 iSNVs (3-106 per sample) were identified and distributed in a highly individualized pattern. The number of iSNVs paced with infection duration peaked within days and declined thereafter. These iSNVs did not stochastically arise due to a strong bias toward C > U/G > A and U > C/A > G substitutions in reciprocal proportion with escalating disease severity. Eight nonsynonymous iSNVs in the receptor-binding domain could escape from neutralization, and eighteen iSNVs were significantly associated with specific HLA alleles. The level and pattern of iSNVs reflect the in vivo viral-host interaction and the disease pathogenesis.

In Depth Viral Diversity Analysis in Atypical Neurological and Neonatal Chikungunya Infections in Rio de Janeiro, Brazil.

Chikungunya virus (CHIKV) is an arthropod-borne virus (arbovirus) transmitted by Aedes mosquitoes. The human infection usually manifests as a febrile and incapacitating arthritogenic illness, self-limiting and non-lethal. However, since 2013, CHIKV spreading through the tropics and to the Americas was accompanied by an increasing number of cases of atypical disease presentation, namely severe neuropathies and neonatal infection due to intrapartum vertical transmission. The pathophysiological mechanisms underlying these conditions have not been fully elucidated. However, arbovirus intrahost genetic diversity is thought to be linked to viral pathogenesis. To determine whether particular viral variants could be somehow associated, we analyzed the intrahost genetic diversity of CHIKV in three infected patients with neurological manifestations and three mothers infected during the intrapartum period, as well as their babies following vertical transmission. No statistically supported differences were observed for the genetic variability (nucleotide substitutions/gene length) along the genome between the groups. However, the newborn and cerebrospinal fluid samples (corresponding to virus passed through the placenta and/or the blood-brain barrier (BBB)) presented a different composition of their intrahost mutant ensembles compared to maternal or patient serum samples, even when concurrent. This finding could be consistent with the unidirectional virus transmission through these barriers, and the effect of selective bottlenecks during the transmission event. In addition, a higher proportion of defective variants (insertions/deletions and stop codons) was detected in the CSF and maternal samples and those were mainly distributed within the viral non-structural genes. Since defective viral genomes in RNA viruses are known to contribute to the outcome of acute viral infections and influence disease severity, their role in these atypical cases should be further investigated. Finally, with the in silico approach adopted, we detected no relevant non-conservative mutational pattern that could provide any hint of the pathophysiological mechanisms underlying these atypical cases. The present analysis represents a unique contribution to our understanding of the transmission events in these cases and generates hypotheses regarding underlying mechanisms, that can be explored further.

Intrahost SARS-CoV-2 k-mer Identification Method (iSKIM) for Rapid Detection of Mutations of Concern Reveals Emergence of Global Mutation Patterns.

Despite unprecedented global sequencing and surveillance of SARS-CoV-2, timely identification of the emergence and spread of novel variants of concern (VoCs) remains a challenge. Several million raw genome sequencing runs are now publicly available. We sought to survey these datasets for intrahost variation to study emerging mutations of concern. We developed iSKIM ("intrahost SARS-CoV-2 k-mer identification method") to relatively quickly and efficiently screen the many SARS-CoV-2 datasets to identify intrahost mutations belonging to lineages of concern. Certain mutations surged in frequency as intrahost minor variants just prior to, or while lineages of concern arose. The Spike N501Y change common to several VoCs was found as a minor variant in 834 samples as early as October 2020. This coincides with the timing of the first detected samples with this mutation in the Alpha/B.1.1.7 and Beta/B.1.351 lineages. Using iSKIM, we also found that Spike L452R was detected as an intrahost minor variant as early as September 2020, prior to the observed rise of the Epsilon/B.1.429/B.1.427 lineages in late 2020. iSKIM rapidly screens for mutations of interest in raw data, prior to genome assembly, and can be used to detect increases in intrahost variants, potentially providing an early indication of novel variant spread.

Remdesivir-induced emergence of SARS-CoV2 variants in patients with prolonged infection.

We here investigate the impact of antiviral treatments such as remdesivir on intra-host genomic diversity and emergence of SARS-CoV2 variants in patients with a prolonged course of infection. Sequencing and variant analysis performed in 112 longitudinal respiratory samples from 14 SARS-CoV2-infected patients with severe disease progression show that major frequency variants do not generally arise during prolonged infection. However, remdesivir treatment can increase intra-host genomic diversity and result in the emergence of novel major variant species harboring fixed mutations. This is particularly evident in a patient with B cell depletion who rapidly developed mutations in the RNA-dependent RNA polymerase gene following remdesivir treatment. Remdesivir treatment-associated emergence of novel variants is of great interest in light of current treatment guidelines for hospitalized patients suffering from severe SARS-CoV2 disease, as well as the potential use of remdesivir to preventively treat non-hospitalized patients at high risk for severe disease progression.

SARS-CoV-2 intra-host evolution during prolonged infection in an immunocompromised patient.

OBJECTIVES: Intra-host SARS-CoV-2 evolution during chronic infection in immunocompromised hosts has been suggested as being the possible trigger of the emergence of new variants. METHODS: Using a deep sequencing approach, we investigated the SARS-CoV-2 intra-host genetic evolution in a patient with HIV over a period of 109 days. RESULTS: Sequencing of nasopharyngeal swabs at three time points demonstrated dynamic changes in the viral population, with the emergence of 26 amino acid mutations and two deletions, 57% of them in the Spike protein. Such a combination of mutations has never been observed in other SARS-CoV-2 lineages detected so far. CONCLUSION: Our data confirm that persistent infection in certain immunocompromised individuals for a long time may favor the dangerous emergence of new SARS-CoV-2 variants with immune evasion properties.

SARS-CoV-2's high rate of genetic mutation under immune selective pressure: from oropharyngeal B.1.1.7 to intrapulmonary B.1.533 in a vaccinated patient.

This is the case report of an 84-year-old man affected by COVID-19 between the 2 doses of vaccination, with negative exitus. We analyzed nasopharyngeal samples of viral RNA collected during the disease and nasopharyngeal and lung samples collected postmortem by reverse transcription LAMP (RT-LAMP) PCR and Next Generation Sequencing (NGS). NGS results were analyzed with different bioinformatic tools to define virus lineages and the related single-nucleotide polymorphisms (SNPs). Both lung and nasopharyngeal samples tested positive for SARS-CoV-2 on RT-LAMP. Through bioinformatic analysis, 2 viral RNAs from the nasal swabs, which belonged to the B.1.1.7 lineage, and 1 viral RNA from the lung sample, which belonged to the B.1.533 lineage, were identified. This genetic observation suggested that SARS-CoV-2 tends to change under selective pressure. The high mutation rate of ORFa1b, containing a replicase gene, was a biological image of a complex viral survival system.

Large-scale analysis of SARS-CoV-2 synonymous mutations reveals the adaptation to the human codon usage during the virus evolution.

Many large national and transnational studies have been dedicated to the analysis of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) genome, most of which focused on missense and nonsense mutations. However, approximately 30 per cent of the SARS-CoV-2 variants are synonymous, therefore changing the target codon without affecting the corresponding protein sequence. By performing a large-scale analysis of sequencing data generated from almost 400,000 SARS-CoV-2 samples, we show that silent mutations increasing the similarity of viral codons to the human ones tend to fixate in the viral genome overtime. This indicates that SARS-CoV-2 codon usage is adapting to the human host, likely improving its effectiveness in using the human aminoacyl-tRNA set through the accumulation of deceitfully neutral silent mutations. One-Sentence Summary. Synonymous SARS-CoV-2 mutations related to the activity of different mutational processes may positively impact viral evolution by increasing its adaptation to the human codon usage.

Limited genomic reconstruction of SARS-CoV-2 transmission history within local epidemiological clusters.

A detailed understanding of how and when severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission occurs is crucial for designing effective prevention measures. Other than contact tracing, genome sequencing provides information to help infer who infected whom. However, the effectiveness of the genomic approach in this context depends on both (high enough) mutation and (low enough) transmission rates. Today, the level of resolution that we can obtain when describing SARS-CoV-2 outbreaks using just genomic information alone remains unclear. In order to answer this question, we sequenced forty-nine SARS-CoV-2 patient samples from ten local clusters in NW Spain for which partial epidemiological information was available and inferred transmission history using genomic variants. Importantly, we obtained high-quality genomic data, sequencing each sample twice and using unique barcodes to exclude cross-sample contamination. Phylogenetic and cluster analyses showed that consensus genomes were generally sufficient to discriminate among independent transmission clusters. However, levels of intrahost variation were low, which prevented in most cases the unambiguous identification of direct transmission events. After filtering out recurrent variants across clusters, the genomic data were generally compatible with the epidemiological information but did not support specific transmission events over possible alternatives. We estimated the effective transmission bottleneck size to be one to two viral particles for sample pairs whose donor-recipient relationship was likely. Our analyses suggest that intrahost genomic variation in SARS-CoV-2 might be generally limited and that homoplasy and recurrent errors complicate identifying shared intrahost variants. Reliable reconstruction of direct SARS-CoV-2 transmission based solely on genomic data seems hindered by a slow mutation rate, potential convergent events, and technical artifacts. Detailed contact tracing seems essential in most cases to study SARS-CoV-2 transmission at high resolution.

Rapid SARS-CoV-2 Intra-Host and Within-Household Emergence of Novel Haplotypes.

In February 2020, the municipality of Vo', a small town near Padua (Italy) was quarantined due to the first coronavirus disease 19 (COVID-19)-related death detected in Italy. To investigate the viral prevalence and clinical features, the entire population was swab tested in two sequential surveys. Here we report the analysis of 87 viral genomes, which revealed that the unique ancestor haplotype introduced in Vo' belongs to lineage B, carrying the mutations G11083T and G26144T. The viral sequences allowed us to investigate the viral evolution while being transmitted within and across households and the effectiveness of the non-pharmaceutical interventions implemented in Vo'. We report, for the first time, evidence that novel viral haplotypes can naturally arise intra-host within an interval as short as two weeks, in approximately 30% of the infected individuals, regardless of symptom severity or immune system deficiencies. Moreover, both phylogenetic and minimum spanning network analyses converge on the hypothesis that the viral sequences evolved from a unique common ancestor haplotype that was carried by an index case. The lockdown extinguished both the viral spread and the emergence of new variants.

Possible recombination between two variants of concern in a COVID-19 patient.

We identified an individual who was coinfected with two SARS-CoV-2 variants of concern, the Beta and Delta variants. The ratio of the relative abundance between the two variants was maintained at 1:9 (Beta:Delta) in 14 days. Furthermore, possible evidence of recombinations in the Orf1ab and Spike genes was found.

Characteristics of SARS-CoV-2 transmission in a medium-sized city with traditional communities during the early COVID-19 epidemic in China.

The nationwide COVID-19 epidemic ended in 2020, a few months after its outbreak in Wuhan, China at the end of 2019. Most COVID-19 cases occurred in Hubei Province, with a few local outbreaks in other provinces of China. A few studies have reported the early SARS-CoV-2 epidemics in several large cities or provinces of China. However, information regarding the early epidemics in small and medium-sized cities, where there are still traditionally large families and community culture is more strongly maintained and thus, transmission profiles may differ, is limited. In this study, we characterized 60 newly sequenced SARS-CoV-2 genomes from Anyang as a representative of small and medium-sized Chinese cities, compared them with more than 400 reference genomes from the early outbreak, and studied the SARS-CoV-2 transmission profiles. Genomic epidemiology revealed multiple SARS-CoV-2 introductions in Anyang and a large-scale expansion of the epidemic because of the large family size. Moreover, our study revealed two transmission patterns in a single outbreak, which were attributed to different social activities. We observed the complete dynamic process of single-nucleotide polymorphism development during community transmission and found that intrahost variant analysis was an effective approach to studying cluster infections. In summary, our study provided new SARS-CoV-2 transmission profiles representative of small and medium-sized Chinese cities as well as information on the evolution of SARS-CoV-2 strains during the early COVID-19 epidemic in China.

Two-step fitness selection for intra-host variations in SARS-CoV-2.

Spontaneous mutations introduce uncertainty into coronavirus disease 2019 (COVID-19) control procedures and vaccine development. Here, we perform a spatiotemporal analysis on intra-host single-nucleotide variants (iSNVs) in 402 clinical samples from 170 affected individuals, which reveals an increase in genetic diversity over time after symptom onset in individuals. Nonsynonymous mutations are overrepresented in the pool of iSNVs but underrepresented at the single-nucleotide polymorphism (SNP) level, suggesting a two-step fitness selection process: a large number of nonsynonymous substitutions are generated in the host (positive selection), and these substitutions tend to be unfixed as SNPs in the population (negative selection). Dynamic iSNV changes in subpopulations with different gender, age, illness severity, and viral shedding time displayed a varied fitness selection process among populations. Our study highlights that iSNVs provide a mutational pool shaping the rapid global evolution of the virus.

Atypical Prolonged Viral Shedding With Intra-Host SARS-CoV-2 Evolution in a Mildly Affected Symptomatic Patient.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection is caused by a respiratory virus with a wide range of manifestations, varying from asymptomatic to fatal cases, with a generally short outcome. However, some individuals present long-term viral shedding. We monitored 38 individuals who were mildly affected by the SARS-CoV-2 infection. Out of the total studied population, three (7.9%) showed atypical events regarding the duration of positivity for viral RNA detection. In one of these atypical cases, a previously HIV-positive male patient presented a SARS-CoV-2 RNA shedding and subgenomic RNA (sgRNA) detected from the upper respiratory tract, respectively, for 232 and 224 days after the onset of the symptoms. The SARS-CoV-2 B.1.1.28 lineage, one of the most prevalent in Brazil in 2020, was identified in this patient in three serial samples. Interestingly, the genomic analyses performed throughout the infectious process showed an increase in the genetic diversity of the B.1.1.28 lineage within the host itself, with viral clearance occurring naturally, without any intervention measures to control the infection. Contrasting widely spread current knowledge, our results indicate that potentially infectious SARS-CoV-2 virus might be shed by much longer periods by some infected patients. This data call attention to better adapted non-pharmacological measures and clinical discharge of patients aiming at preventing the spread of SARS-CoV-2 to the population.